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Abstract 

Integration of nonlinear dynamical systems is usually seen as associated to a symmetry 
reduction, e.g. via momentum map. In Lax integrable systems, as pointed out by 
Kazhdan, Kostant and Sternberg in discussing the Calogero system, one proceeds in 
the opposite way, enlarging the nonlinear system to a system of greater dimension. We 
discuss how this approach is also fruitful in studying non integrable systems, focusing 
on systems in normal form. 

1 Introduction. Dimension increase and splitting 

Finite dimensional integrable systems have been a rarity for a very long time; up to few 
decades ago, they could be easily enumerated: harmonic oscillators, the Kepler system, a 
few spinning tops. 

Apart from harmonic oscillators ~ which give linear equations and are the prototype 
of integrability - these are integrated by exploiting their symmetries and the associated 
integrals of motion to reduce the problem to a lower dimensional one; thus, we have 
symmetry reduction. 

The situation changed radically in the seventies, when entire classes of new integrable 
systems appeared, integrable by the Lax pair construction Among these a prominent 
role - also due to priority - is taken by the Calogero system In this case, one has an 
arbitrary number N of points on the line, interacting via a certain pair potential; thus the 
system is described by a natural hamiltonian H = (1/2) + ^(^)- Calogero showed 

that there is an invertible map from to GL{N), mapping the vector x to a matrix 
L(x), such that the evolution of x{t) results in an evolution of L{t) := L[x{t)] governed by 
Lax equations; thus one integrates L(t) and obtains x{t). 

As pointed out by Kazhdan, Kostant and Sternberg ^Hj, this also initiated a new 
way of integrating systems (which they studied geometrically): indeed, rather than trying 
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to reduce the dimension of the system, Calogero passed to consider a system of higher 
dimension}. In this note we want to show that this approach is fruitful also beyond the 
realm of integrable hamiltonian systems. 

The key observation made by Kazhdan, Kostant and Sternberg ^^1) see also jl5l I16| . 
in the context of Hamiltonian systems is that given a differential equation A in dimension 
n, in some cases it may be helpful to increase its dimension: this is in particular the case 
when it is possible to describe A as originating from either the symmetry reduction or 
the projection to an n-dimensional manifold of a simpler (e.g. linear) equations A* in 
dimension n + m. 

It was remarked in ^01 that this approach does apply also to non-hamiltonian systems, 
and in particular to systems in Poincare-Dulac normal form ^ [31 [20] (a remark to this 
effect was already contained in [201) ■ Essentially, a system in normal form with only a 
finite number of resonances can always be mapped to a linear system, and is equivalent 
to the latter on a certain invariant manifold. As observed there, such a procedure is 
geometrically interesting, but equivalent to a classical technique already known to Dulac 
[S], and abscribed by him to Horn and Lyapounov. In more general terms, i.e. outside 
the scope of normal forms theory, this corresponds to the situation of a finite dimensional 
centralizer for a certain algebra associated to the vector field ,11. 

In the present note, we show how "dimension increasing" nicely combines with a "split- 
ting" (in a sense to be described below) approach [H] in the case of an infinite number of 
resonances^; this case cannot be tackled by the Horn-Lyapounov-Dulac approach. 

Roughly speaking, the procedure presented in this note combines the two approaches 
mentioned above, i.e. describing a nonlinear system as projection of a linear one and 
symmetry reduction in the sense described in [Hj. 

A more abstract (and general) treatment is given elsewhere |12j : in the present note 
we adopt an approach and notation aimed at applications, and discuss a number of con- 
crete applications. We focus on systems in normal form (sect. 2) and in particular on the 
embedding and splitting of systems with infinitely many resonances (sect. 3). Some simple 
examples are discussed in sect. 4, while sect. 5 is devoted to (in general non-hamiltonian) 
nonlinearly perturbed oscillators, and sect. 6 to bifurcation problems. 



2 Equations in normal form; resonances 

We consider an ODE in i?" with a fixed point in the origin and expanded around this in 
a power series; we write this in the form 

J oo 

^ = fix) = Ax + ^Mx) (1) 

k=l 

where is polynomial with fk{ax) = a^'^^^x), and we have singled out the linear part 
/o(x) = Ax. We also write F for the nonlinear part of /, i.e. /*(x) = AjX^ + F^{x). 

^Another classical instance of simplification by dimension increase is the trick to linearize the matrix 
Riccati equation. See also |11 118| in this context. 

^For a different use of the splitting approach in normal forms type problems, see \'22\. 
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As well known, the matrix A can be uniquely decomposed into its semisimple and 
nilpotent parts, which commute with each other and hence with A: 

A = Ag + An ; [As, An] = 0. 

In the following, we will use the following vector fields {di := d/dx^): 

Xa := (AxYd^ , Xo = {Asfd^ , Xp = F\x)d^ , Xf = f{x)d^ . 



2.1 Resonances 

Let us denote by {Ai, A„} the eigenvalues of A; take a basis {ei, e„} in i?" consisting 
of generalized eigenvectors of A, i.e. eigenvectors of Ag: A^ej = Xjej. We will use x 
coordinates in this basis, and the multiindex notation 

We say that the vector monomial v^^^ := x^e^ is resonant with A if 

n n 

(/U • A) :=^//iAi = Xa with fii > , |^|:=^/ii>l. (2) 

i=l i=l 

The relation • A) = A^ is said to be a resonance relation related to the eigenvalue Xa, 
and the integer \n\ is said to be the order of the resonance. In our context it is useful to 
include order one resonances in the definition (albeit the trivial order one resonances given 
by Xa = Xa are of little interest). Note that here one could as well consider Ag rather 
than A. 

The space of vectors resonant with (the semisimple part of) A is defined as the linear 
span of the vectors v^^q, defined above. 



2.2 Normal forms 

We say (see e.g. ^EH^O]) that (1) is in Poincare-Dulac normal form^ if its nonlinear part 
F{x) is resonant with A. This implies that 

[Xo,Xf] = . (3) 

It should be mentioned that the presence of a nilpotent part An in A introduces some 
subtleties. If (3) holds, then both Xa and Xp commute with Xq, and therefore 

[^0,^/] = 0, (4) 

i.e. the system has a symmetry described by a semisimple matrix. 

As well known, starting from any dynamical system (or vector field) of the form (1), we 
can arrive at a dynamical system (or vector field) in Poincare-Dulac normal form by means 
of a sequence (in general, infinite) of near-identity transformations obtained by means of 
the Poincare algorithm; these combine into a near-identity transformation H defined by a 
series which is in general only formal. 

^The reader should be warned that some different definition is also in use 0. 



4 



G Gaeta, S Walcher 



Remark 1. We may reformulate the definition of systems in normal form by saying that 
/ is is normal form if and only li Xp is in the centralizer of Xq. 

Remark 2. If the system is required to have some symmetry, say [Xf,Xg] = with (in 
an obvious notation) g^{x) = BjX^ + G*, then = as well, i.e. Xp is in the 

centralizer s of both Xq and Xb- More generally, if ^/ is an element of some Lie algebra 
Q, then under suitable conditions, see 0, it can be put in joint normal form, i.e. (again 
with obvious notations) [Xj^Xb^] = 0. 

2.3 Sporadic resonances and invariance relations 

Let us consider again the resonance equation (2). It is clear that if there are non- negative 
integers ai (some of them nonzero) such that 

n 

^ CT^\^ = , (5) 
i=l 

then we always have infinitely many resonances. The monomial cj) = will be called a 
resonant scalar monomial. It is an invariant of Xq, and any multiindex fi with fii = kai+5ia 
provides a resonance relation {fJ. ■ X) = Xa related to the eigenvalue A^; in other words, 
any monomial x'^^'x" = (/>x" is resonant, and so is any vector ^ki7-^ea,(y' 

Therefore, we say that (5) identifies a invariance relation. The presence of invariance 
relations is the only way to have infinitely many resonances in a finite dimensional system 
(see 1201). 

Any nontrivial resonance (2) such that there is no a with o"j < fii (for all i = 1, ...,n) 
providing an invariance relation, is said to be a sporadic resonance. Sporadic resonances 
are always in finite number (if any) in a finite dimensional system j2()j . 

Any invariance relation (5) such that there is no i' with < Cj (and of course v ^ 
providing another invariance relation, is said to be an elementary invariance relation. 
Every invariance relation is a linear combination (with nonnegative integer coefficients) of 
elementary ones. Elementary invariance relations are always in finite number (if any) in 
a finite dimensional system j2()j . 

3 Embedding systems with invariance relations in quasi- 
linear systems 

In this section we will discuss how the procedure described in JOl generalizes, in connection 
with a "splitting" of the system described by y [S], in the presence of invariance relations. 

We should preliminarily identify all sporadic resonances (^ • A) = A^ and elementary 
invariance relations (o" • A) = 0. We associate resonant monomials x^ and resonant vectors 
^pi.,a to the former ones, and invariant monomials x'^ to the latter ones. 

We then introduce two set of new coordinates: these will be the coordinates tf^, tu*" in 
correspondence with sporadic resonances (as in |10j^. and other new coordinates 0^, 
in correspondence with elementary invariance relations. 
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We should also assign evolution equations for the w and coordinates; these will be 
given in agreement with (1) itself. That is, the equations for the w will be 

div^ dnj"^ dx^ 
and as for the 0's we assign 

d</>" dcj)"" dx* 



z''ix,w,ct>) . (7) 



We will thus consider the enlarged space W = {x,w,(j)) = and in this the 

vector field 

Y = f{x,w,4>)^, + h\x,wA)^ + z^{x,w,4>)^^. (8) 

Note that some ambiguity is present here, in that we can write the coefficients of this 
vector field in different ways as a function of the x, w, (p. Indeed, the vector field Y is 
uniquely defined only on the manifold identified by ^/^* := w"^ — x'^(') = 0, (j)a — Ca{x) = 0. 

Lemma 1. The {n + m)- dimensional manifold M C W identified by ij)^ := — x^ = 
is invariant under the flow ofY. 

Proof. Obvious by construction. A 

Lemma 2. The functions defined in (7) can be written in terms of the (j) variables 
alone, i.e. dz"- /dx^ = dz"" jdw^ = 0. 

Proof. Every analytic invariant of Xq can be represented as a convergent series in the 
(see |2n|)i their evolution is also invariant under the Xq action, hence can also be written 
in terms of invariants, hence of the z"" themselves. A 

Corollary 1. The evolution of the (j) variables is described by a (nonlinear) equation in 
the (j) variables only. 

Proof. This is merely a restatement of lemma 2 above. Note that the equations for x 
and w depend on (p and are therefore nonautonomous. A 

Proposition 3. The analytic functions /* and h^ defined above can be written as linear 
in the x and w variables, the coefficients being functions of the (p variables. 

Proof. Recall each is a monomial = x^ ...x^ with [p, A) = for some s = 1, n; 

is) 

with reference to this integer s, we add a label to Wj, i.e. write w- . By construction and 

by the results above, each /"^ can be written in the form /"^ = a™((/)) -Xm + X^fc c™((/>)t(;^™'^ , 
with analytic am and c™. So the assertion for the is obvious. 

The time derivative of ti'j*'' under the flow of (1) will be wj*^ = {dwj^^/dx"^)f"^. There- 
fore it is sufficient to show that (dwj^^ / dx"^)w^j!"^ is zero or a multiple of x™" or of some vu^'^^ 
with a suitable resonant scalar monomial as a factor. But the above operation just means 
to replace one factor x™ by w^^^ , and the resulting linear combination of the eigenvalues 
still yields a resonance relation (2) with A^ on the right hand side. See jl2j for a different 
approach to the proof. A 
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Corollary 2. The evolution of the x and w variables is described by nonautonomous 
linear equations, obtained by inserting the solution 4> = 4>{t) of the equations for cj) in the 
general equations x = f{x,w,(l)), w = h{x,w,(l)). 

Proof. Obvious. A 

Remark 3. We will also say that the vector field Y is quasi-linear, meaning by this that 
it is linear in the x and w variables. In this way we recover - as a special case - the 
situation discussed in M as well as the terminology used there. 

Remark 4. The results obtained here extend and unify those given in [HIE]; see also |2U| . 
As the (j) identify group orbits for the group G generated by the Lie algebra, we interpret 
(j) = z{(l)) as an equation in orbit space, and the equation for (x, w) as an equation on the 
Lie group G. Methods for the solution of the latter are discussed in |2J, see also j^- 

Remark 5. If no invariance relations are present, hence no (j) variables are introduced, 
then the system describing the time evolution of the x, w variables is linear; this is the 
situation studied in [Tn|. Note that in this case we have exactly the interpretation of 
normal forms as projection of a linear system to an invariant manifold, without symmetry 
reduction. 

Remark 6. If there are no sporadic resonances of order greater than one then Proposition 

3 yields a linear system for the /*, with functions of the (p variables as coefficients. There- 
fore, upon solving the reduced equation for the (p variables one obtains a non-autonomous 
linear system. Moreover, if all eigenvalues are distinct then we have a product system of 
one-dimensional equations. 

Remark 7. Finally, we note that if (f){t) converges to some (f)o (this is always the case 
if the (j) space is one-dimensional and \4>it)\ does not escape to infinity), the asymptotic 
evolution of the system is governed by a linear autonomous equation for x and w (see ^Hl 
for the behavior of asymptotically autonomous equations). Similarly, if there is a periodic 
solution (f)(t) with (p{t) — s- (f){t), the asymptotic evolution of the system is governed by a 
linear equation with periodic coefficients for x and w. 

4 Examples 

Example 1. (See ^Ol)- For A = diag(l, k), A; G N, the only resonant vector is v = x^e2, 
corresponding to a sporadic resonance, and there is no invariance relation. Systems in 
normal form correspond to 

X = xdx + {ky + cx'^)dy 
with c a real constant. According to our procedure, we define w = x^, and obtain 

Y = xdx + (ky + cw)dy + kwd^ ', 

the invariant manifold M is given hy Tp := w — x'' = 0. The solution to the system in W for 
initial data {xo,yo,WQ) is x{t) = xqc*, y{t) = yoe^^ + {cikwQ)te^* , w{t) = WQe^^; for initial 
data on M, i.e. wq = Xq, the solution remains on M and its projection to R^ = {x,y) is 
x{t) = xoe*, y{t) = [yo + (ci/cxg)t]e'=*. 
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Example 2. Consider the matrix 




with eigenvalues 1). There is one elementary invariance relation, Ai + A2 = 0, and 

no sporadic resonance. The linear centralizer is spanned by A itself and by matrices Di = 
diag(l, 1, 0) and D2 = diag(0, 0, 1). The ring of invariants is generated by := + y^. 
Systems in normal forms are written as 

X = a(r^)x — P{r'^)y 
y = (5{r'^)x + a{r^)y 
z = ^{r'^)z 

where a, /?, 7 are arbitrary power series. 

Following our procedure, we introduce one further variable cp = r"^; for this we have 
d(f)/dt = 2{xx + yy) = 2r^a(r^). Hence the system in = R"^ is written as 

X = a{(j))x — P{(j))y 

y = (3{(t))x + a{(j))y 
z = j{(t))z 
(j) = 2(t)a{(t)) . 



Example 3. Consider the matrix 



A 



1 



Vo 










0\ 



1 

/c/ 



(A; G N, A; > 1) with eigenvalues (— +?', 1, k). There is one sporadic resonance, = A4, 
and one elementary invariance relation, Ai + A2 = 0. The linear centralizer is spanned the 
matrices Di = diag(l, 1, 0, 0), D2 = diag(0, 0, 1, 0), D3 = diag(0, 0, 0, 1) together with 



Ml 



/O -1 0\ 

10 



yo 00/ 



and M2 = 



1^ 


































Vo 





k 





Systems in normal forms are written as 

xi = a(r^)xi — P{r'^)x2 
X2 = /?(r^)xi + a(r^)x2 
X3 = 'yir'^)x3 
±4 = r]{r'^)x4 + 9{r'^)x2 



where := + y^ and a, j3, 7, 77, 6 are arbitrary power series. 
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Following our procedure we introduce = r^, for which dcp/dt = 2{xx+yy) = 2r^a(r^), 
and w = for which dw/dt = kx'^~^X3 = kj{r'^)x^. Hence the system in = is 
written as 

xi = a{(t>)xi - P{(t>)x2 
X2 = /3{(f>)xi + a{4>)x2 

w = k^{(t>)w 

5 Perturbed oscillators 

Example 4. Perturbation of oscillators in 1:1 resonance. Consider the matrix 

/O -1 \ 
10 
-1 ' 

Vo 1 / 

which we also write in block form as 

with eigenvalues There are no sporadic resonances of order greater than 

one, and four elementary invariance relations: 

Ai + A2 = , A3 + A4 = , Ai + A4 = , A2 + A3 = ; 

all other resonances can be described in terms of these. We stress that the equations 
describing these invariance relations are linearly dependent; however they should be con- 
sidered, according to our definition, as different elementary ones. Corresponding to this, 
the associated invariant quantities will not be functionally independent (obviously, we 
cannot have more than three independent invariants for a flow in R^). 

The linear centralizer of A is an eight-dimensional algebra, spanned by the following 
matrices (in two by two block notation): 

^^ = (0 0) ' ^2 = (o /) ' ^^=(1 0) ' ^'={-J 0) ' 

(9) 

^'-[0 oj ' Vo j) ' ^'-[-i oj ' Vj 0) ■ 

Note here we have chosen a basis with Bi = Bf, Si = —S^. 

The linear system ^ = describes two oscillators in 1:1 resonance; the normal form 
will correspond to a perturbation of these, generically breaking the exchange symmetry 
among the two oscillators. Systems in normal form are compactly written as 





/ a 


-/5 


7 


-v\ 


/x\ 


y 




a 


V 


1 


y 


z 


/i 


—v 


a 


— T 


z 


\wj 




1^ 


T 


^ J 


\wj 
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where a,/3, ...,r are arbitrary power series in the elementary invariants 

4)1 = + , 4>2 = + w'^ , (f)3 = xz + yw , 4>4 = xw - yz ■ 

note that (f)a = {^,Ba^), with (., .) the scalar product. We abbreviate the above evolution 
equation as 

i = m) i ■ 

The evohition equations for the as required by our procedure, are simply (only the 
selfadjoint matrices defined in (9) appear) 

Example 5. Perturbation of oscillators in 1 : A; resonance. Consider the matrix 



A = 



/O -1 \ 

10 

-A; 

Vo fc y 



(with /c S N, /c > 1) with eigenvalues {—i,+i,—ik,ik). This is put in diagonal form 
passing to variables 

Ci = {xi - ix2)/2 , ^2 = {x\ + ix'2)/2 , ^3 = (xs - 1x^)12 , ^4 = (X3 + ix^j/i , 
which we use in intermediate computations below. We also write ^ = Ax, x = A""-*^^, with 



A = 



/I -i \ 

110 

1 -i 

Vo 1 i y 



; A 



-1 



/I 1 \ 

i -i 

11 

Vo i -ij 



There are four sporadic resonances: 

kXi = A3 , kX2 = A4 , A3 + (fc — 1)A2 = Ai , A4 + (fe — l)Ai = A2 . 

Moreover, there are four elementary invar iance relations: 

Ai + A2 = , A3 + A4 = , AAi + A4 = , fcA2 + A3 = . 

All other resonances can be described in terms of these. 
Systems in normal forms are written as 

ii = ai6 + Mt^Cs 

6 = "26 + Mi'^U 

U = 0:4^4 + '&4^2 > 



(10) 



where Q;i,t?i are arbitrary power series in the invariants of the linear fiow. 
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Following our procedure, we introduce variables 

related to sporadic resonances. We also introduce variables related to elementary invari- 
ance relations, given by 

01 = 66 , <t>2 = 66 , h = CiU , (f>4 = C26 • 

We must then introduce evolution equations for the w and (f> variables according to our 
procedure, i.e. according to (6) and (7) above. As for the w, wc get 

Wl = kaiwi + i)i4'i^^(,3 

W2 = ka2W2 + ^24>i~^U 

W3 = [as + {k- I)a2]w3 + [{k - l)-&2(t>i~^ h + ^30i"^]6 
W4 = [a4 + ik- I)a2]w4 + [{k - l)^i^\-'^^2 + 1?4'^^^]6 • 

Let us now consider the equations for the 0; we easily get 

4>i = (ai + a2)<Pi + ^?203 + '&i4>A 

^3 = (kai + ai)(l)3 + k^i(l)\-^<l)2 + '&4<lA 
^4 = {ka2 + q;3)'?!'4 + A;i?2'?!'i~ V2 + '&3<t>i ■ 

Summarizing, all systems of the form (10) - i.e. in normal form with respect to the 
linear part ^ = — are written as the autonomous system (11) for the (j) variables, plus 
a linear nonautonomous system, which introducing the notation rj = {i;w) can be written 

as 

?7 = M 77 

where M = M^cf)) is a matrix which we write as M = D + L, where 

D = di&g{ai,a2,a3,a4;kai,ka2,a3 + {k - I)a2,a4 + {k - l)ai) 
and L is an off-diagonal sparse matrix with nonzero terms 

Li7 = , L28 = ^2 , L35 = 1)3 , L46 = 1?4 ; L53 = , L64 = ^2(l)i~^ , 

Ln = [{k - l)^?20^'</'2 + ^34>i~^] , Ls2 = [{k - l)^i4-''4>2 + ^a4>i~^] ■ 

Example 6. Perturbation of two oscillators with no resonance. Consider the 
matrix 

/O -1 \ 
10 
-TT 

yo TT / 

with eigenvalues {—i, +i, —iri, m). 

Now there are no sporadic resonances of order greater than one, and two elementary 
invariance relations: 

Ai + A2 = , A3 + A4 = . 
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Note that now we have only two invariants, rl = x\+ X2 and r| = x| + x|: this is easily 
understood as we have irrational flow on the two-torus C R^. The centralizer of A 
corresponds to matrices 

/a -b \ 
6 a 
c -d 
\(} Q d c J 

Thus systems in normal form will be written as 



M 







/ a 


-P 





^ 






±2 




(3 


a 












X3 










7 








\X4) 




^0 





V 


7 J 




\Xil 



with a, 13, 7, ?7 being power series in , r|. 

Following our procedure we introduce variables 0i = + a;| and 
evolution is given by 

^1 = 2q;^i , ^2 = 27^2 • 



5>2 = x\ + x\, whose 



Example 7. Perturbation of oscillators in 1:1:1 resonance. Consider the six- 
dimensional matrix written in block form as 



A = 



to J 












uJ 




where J 








uJ 





-1 



Passing to coordinates {CjyVj) = {Pj + Hj^Pj ~ tliis reads 

A = dia,g{iLO, —iLO,iu!, —iuj,iLO, —iu) . 

The eigenvalues are Aj^ = {—l)''iuj, hence there is no sporadic resonance of order greater 
than one^ and nine invar iance relations: 

Ai + A2 = A3 + A4 = A5 + Ae = 
Ai + A4 = Ai + Ae = A3 + Ag = 
A2 + A3 = A2 + A5 = A4 + A5 = 

The invariants corresponding to these are of course 

^1 = ^im, h = 6%, ^3 = 6%, 

Going back to the original coordinates, expressions for these are recovered from (no sum 
on repeated indices) 

C,Vj = {p] + ei) ; 0% = {p,Pk + qm)+^{q,Pk-p,qk) 
''Note these are present in the case oil : k : I resonance (with integers k,i> 1). 
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and is thus more convenient to pass to a different basis for invariant functions, i.e. 



H = (Pi + Qi) , 'P2 ■■= 4>2 = {pI + ql) 

(pi + ^7)/2 = {piP2 + qiq2) , 4>5 
<k + ^9)/2 = (piP3 + qm) , h 

<i>6 + ?8)/2 = (P2P3 + 9293) , 4"^ 



^, := (pz = {pl + qfj ; 
^4 - ^7)/(2i) = (giP2 -^192) 
^5 - ^9)/(2i) = (giP3 -pm) 
h - ^8)/(2i) = (g2P3 -^293) 



The centrahzer C{A) of A in Mat(6,i?) is an algebra spanned by eighteen matrices; 
these are written in 2 x 2 block form as 

Cii C12 CisX 

C21 C22 ^^23 j 

, C31 C32 C33 / 

where the Cij are two-dimensional matrices written as {aij and (3ij real constants) 



aij I + J 



It is easy to extract from these a basis made of nine selfadjoint matrices B]- = and 
nine antiselfadjoint ones = — S"^- With x = (pi, (?i,P2, 92,P3, 93), and (., .) the standard 
scalar product in these can be chosen so that (j)a = (x, S^x). 

In this compact notation the 1:1:1 resonant three-dimensional oscillator is described 
by X = ^x. We write the normal form of evolution equations for perturbation of this as 



X 



K{4>) X 

(X, {BaK + K+Ba)^) 



with K an arbitrary matrix in C{A). 



6 Bifurcations 

Example 8. Hopf bifurcation. Consider the matrix 

- - C T) 

with eigenvalues Ai = —iooo, A2 = iwo- There is no sporadic resonance, and one elementary 
invariance relation, Ai + A2 = 0, with associated invariant cf) = + . The most general 
system in normal form is therefore x = a{(j))x — f3{(j))y, y = (3{4')x + a{(f))y. According to 
our procedure, the evolution equation for the new coordinate (j) will he (j) = 2(j)a{(j)). As 
the linear part of the system is given by A, we must require a(0) = 0, /3(0) = ujq. 

In applications, one is interested in the case where the system does also depend on an 
external ("control") parameter ^u, which usually does not evolve in time^, the linear part 
being given by A at the critical value. In this case the normal form and the (f) evolution 
equation read 

X = a{(t),n)x - P{cl),i2)y 
y = P{(l),n)x + a{(j),^j.)y 
(j) = 2(/)a(0, /x) . 
different framework is provided by dynamic bifurcations pill7|. 
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In the standard model of Hopf bifurcation, a((p,iJ,) = n — ccp, and we write f3((p,iJ,) = 
ujo + b{<p, jj) with 6(0, 0) = 0. This corresponds to the normal form 

X = fix — uqu — 6(x^ + y^, fi)y — c(x^ + y^)x 
y = ujqx + ny + b{x'^ + y^, - c{x'^ + y'^)x 

which in our approach reads 

x = fix- (Joy - K'P, i^)y - c{(p)x 

y = lvqx + fiy + b{(p, ^x)x - c{(p)x 
^ (j) = 2 0a(0,/i) . 

Note that the space of invariants is one-dimensional (with the additional constraint > 0) ; 
thus, either ^(t) is unbounded for t > 0, or it approaches one of the zeroes of the function 
a{(f),iJ,), say ^o- In this case, the system reduces asymptotically to a linear one: 












The standard analysis of Hopf bifurcation is readily recovered in this way. 

Note that we can also look at Hopf bifurcation in a slightly different way, i.e. include 
the fi variable from the beginning. In this case the matrix A is given by 



A 



with eigenvalues {—i, i, 0) and invariance relations Ai + A2 = and A3 = 0; the associated 
invariants are (p and /x. The linear centralizer of A is spanned by matrices 



M 



and correspondingly the normal form will be 

X = a{(l),^i)x - (3{(l),^)y 
y = I3{(t),ij)x + a{4>,ii)y 

The equation for (f) is just the one given above, and we are led to the same system; 
interpreting jjL as an external control parameter forces ^{(j),iJ,) = (if not, there is a 
feedback of the system on the control parameter). 

Example 9. Hamiltonian Hopf bifurcation. Consider the matrix 




A = 





Vo 



—LO 











\ 


— U! 



fil + uJ 






— /x/ + UjJ 
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with eigenvalues Ai = —n — ico, A2 = —ji + iu), A3 = jJL — ioj, A4 = jJL + iu). We assume w 7^ 
(so it could be rescaled to a; = 1) and n 0; the case n = corresponds to example 4. In 
applications, one considers the case where /x is an external control parameter and studies 
the changes as this is varied; /i = is a critical value. 

The matrix A is diagonalized, for all n, passing to variables as in example 5 above. 
The evolution x = Ax preserves the symplectic structure k = dx^ A dx^ + dx^ A dx^. 

It is easy to check, for generic that there is no sporadic resonance and that there 
are two elementary invariance relations, given by Ai + A4 = and A2 + A3 = 0. The 
corresponding invariants will be (pi = <y92 = they are complex conjugate, 

and correspond to real invariants (f)i = X1X3 + X2X4 and (j)2 = X1X4 — X2Xs- These are 
also written for later reference as (f)a = (l/2)(^, 5^^), with Ba obvious four-dimensional 
symmetric matrices. 

The linear centralizer of A (for 7^ 0) is given by matrices written in block form, with 
ak and Pk real constants, as 



where now ak,(3k will be functions of the invariants (j)i,4>2- Note that such systems in 
general do not preserve the symplectic form k, unless a2 = — ai and (32 = Pi- The 
evolution of the (p is given by 4>a = (1/2) {x, {M^Ba + BaM)x). 

It is convenient to write the system in terms of the two-dimensional vectors r]i = 
{x^,x^), r]2 = {x^,x^) and (/) = {(/)i,(l)2)- Moreover, we single out the linear part in 
the functions aj. and Pk, writing ak{(t>) = {—l)^^^fJ, + a{(f)), Pk{4') = ^ + bi.{(f)) with 
'ifc(O) = ^A;(0) = 0; the smooth functions are a^, bk are arbitrary apart from this constraint, 
and can also depend on the parameters /x and u. 

The system is hence described in our approach by the following equations: 



If ak,bk are such that the system is hamiltonian, the (f) are always strictly invariant, 
and we are reduced to a linear system on each level set of (f); if the a^jb^ are such that 
the system is not hamiltonian but the (two-dimensional) equation for the ^ satisfies the 
condition for the existence of a limit cycle, as is often the case in bifurcation problems, 
then remark 7 applies. 

Let us look more closely to the perturbation of the case /x = 0; as recalled above, the 
analysis of this case can be recovered from example 4. More precisely, we can rescale time 
so that (J = 1, and allow the arbitrary functions appearing in the analysis of example 4 to 
also depend on the parameter /i. 

Alternatively, we can proceed as at the end of the previous example, and include the /Li 
variable from the beginning. With a 2 © 2 1 block notation, we have now 




Correspondingly, systems in normal form will be given by 



X = Mx 



771= {nl + ujj)r]i + [oi ((/))• I - 61 ((/))• J] r?i , 

m = {-fj-i + ^J) m + [a2i(l>) ■ I - b2{4>) ■ J] m , 

^ = [{aiicP) + 02(0)) • / + M^) - 6i(0)) • J] . 



A = 





-fil + ujJ 
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with linear centralizer 

'aiI + l5iJ 0^ 

M = I a2/ + /32J 


In the normal form, going back to the original time parametrization, we have 

??i = (iJ.1 + ioJ)r]i + [ai{(j),iJ:) ■ I - f3i{(p,iJ,) ■ J]r]i , 
m= {-nl + uJ)r]2 + [q;2(^,/x) • / - /?2(0,Ai) • <^]»?2 , 

Again, interpreting ji as an external control parameter requires 7((?!), ji) = 0. 
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